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Abstract 

The pressure of QCD admits at high temperatures a factorization into purely perturba- 
tive contributions from "hard" thermal momenta, and slowly convergent as well as non- 
perturbative contributions from "soft" thermal momenta. The latter can be related to various 
effective gluon condensates in a dimensionally reduced effective field theory, and measured 
there through lattice simulations. Practical measurements of one of the relevant condensates 
have suffered, however, from difficulties in extrapolating convincingly to the continuum limit. 
In order to gain insight on this problem, we employ Numerical Stochastic Perturbation The- 
ory to estimate the problematic condensate up to 4-loop order in lattice perturbation theory. 
Our results seem to confirm the presence of "large" discretization effects, going like aln(l/a), 
where a is the lattice spacing. For definite conclusions, however, it would be helpful to repeat 
the corresponding part of our study with standard lattice perturbation theory techniques. 
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1. Introduction 



Given possible applications in cosmology and in the phenomenology of heavy ion collision 
experiments, as well as the important theoretical role that the free energy density plays in 
understanding the properties of any finite-temperature system, the pressure (or minus the 
free energy density) of Quantum Chromodynamics (QCD) is one of the central observables of 
relativistic thermal field theory. In this paper we focus on its determination at temperatures 
above a few hundred MeV, where the system is deconfined; there, if anywhere, it should 
eventually be possible to establish a quantitative first-principles description in terms of the 
temperature and the fundamental parameters of the theory. 

Given the importance of the problem, a huge amount of work has already been carried 
out on the topic. (The references most directly related to the present work are cited at 
the beginning of the next section.) In general, the approaches can be divided into numerical 
(i.e. lattice Monte Carlo) and analytic (i.e. weak-coupling expansion or various improvements 
thereof) techniques. In addition, there is a strategy — the one that we follow here — which 
combines elements from both sides. The idea is to factorise the system into two parts: "hard" 
momenta, whose contribution is perturbative, and "soft" momenta, which need to be treated 
non-perturbatively. (One benefit of this approach is that dynamical quarks remain cheap 
even in the chiral continuum limit, since only gauge fields possess soft momenta.) Our study 
concerns the non-perturbative soft part, but not directly its lattice measurement; rather, the 
point is that for such a factorization to work, both sides of the result need to be converted 
to the same regularization scheme, so that they can be added together. In some situations, 
like in ours, such a scheme conversion can turn out to be technically as demanding as the 
non-perturbative measurement itself, and this is the ultimate challenge that we try to tackle 
in this work. 

This paper is organised as follows. The general setup of factorising the pressure to contribu- 
tions from perturbatively computable terms, and ones that need to be estimated numerically, 
as well as the role that the present study plays in this setup, is outlined in Sec. [H The 
outline is made quantitative in Sec. [3l where we define the precise quantities that we want to 
determine. The tool used for the computation. Numerical Stochastic Perturbation Theory, 
is reviewed in Sec. [H The numerical data is analysed in Sec. O our results and conclusions 
comprise Sec. [6j In three appendices, we detail the perturbative expressions that we have 
worked out explicitly in lattice regularization. 

2. Outline of setup 

At a high temperature T and a small gauge coupling constant g, there are parametrically 
three different momentum scales in hot QCD: k ~ -kT , gT , g^T / a All the effects of the 
hard scale, k ~ vrT, can be accounted for by a method called dimensional reduction [Tl[2]. In 
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particular, the pressure, or minus the free energy density, can be written as [3] 



Pqcd{T) = pe{T)+ hm -In / PA^PAg expf-^e) 



(2.1) 



where V is the volume, A" are gauge fields, Aq are scalar fields in the adjoint representation, 
and Se is a three-dimensional effective field theory, to be specified presently (the subscript 
may refer to "Electrostatic QCD" [3j). The function PEiT) gets contributions only from 
the hard scale, k ~ vrT, and is computable in perturbation theory; the path integral with 
Se contains the contributions of the soft modes, k ~ gT, g'^T/'K, and should preferably be 
determined non-perturbatively (the contributions from the modes k ~ g^T/ir are genuinely 
non-perturbative [H [5]; those from the modes k ~ gT are in principle still perturbative, 
but in general slowly convergent [3], [6]-[9], although some observables with possibly faster 
convergence have also been found [lO]). The same description applies also in the presence of 
a small quark chemical potential [Tl], allowing to compute further quantities such as quark 
number susceptibilities [l2]. At least when treated non-perturbatively, Electrostatic QCD 
appears to yield a quantitative description of the full four-dimensional theory up from about 
T ~ l.STc, where Tc denotes the pseudocritical temperature of the QCD crossover (see, e.g., 
ref. [13] and references therein). 

Now, to give a precise meaning to Eq. ()2.ip . requires the specification of a regulariza- 
tion scheme. Though a purely perturbative challenge, the determination of Pe{T) is fairly 
complicated in practice [6]. Therefore, it is preferable to use dimensional regularization for 
the algebra: all the 3-loop and 4-loop results available today have been obtained in the MS 
scheme [3], [6]-[8], [Hill]. 

Let us, correspondingly, denote the MS scheme "vacuum energy density" of the theory 
defined by Se, with 



Though each part is scheme-dependent, the expression as a whole is not. In the following, 
we concentrate exclusively on the determination of f ^s- 

To take further steps, we need to specify the effective action Se- It reads 



[V^ooV J V /J_ 

where V = f d'^x is the d-dimensional volume. Then Eq. (|2.ip becomes 




(2.2) 




(2.3) 




(2.4) 



Ce 



- Tr [4] + Tr [A, ^o]' + m^TY [AI\ + A3(Tr [Al]f + ... . 



(2.5) 
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Here i = d = 3 - 2e, Fi^ = {i/ g^)[D,, Dj], Di = 8, - ig^A,, Ai = ^f^^ ^lo = 

AI^T", and are hermitean generators of SU(A^c), normalised as Tt[T'^T^] = 6'^^ /2U The 
dimensionalities of and A3 are GeV^"''^'^. There are also higher order operators, classified 
in ref. [17], whose parametric importance has been analysed in ref. [8]; even though some 
of them contribute at the same parametric order as some of the effects that we are after, it 
is still a well-defined problem to start by determining the full non-perturbative effect of the 
truncated form of the theory in Eq. ()2.5p . Therefore we will ignore all higher order operators 
in the following. 

Now, f^s may include a part which is independent of rrig. However, this part can be 
evaluated by sending mg — > 00, whereby the field Aq can be integrated out. Thereby the 
problem reduces to that already considered in refs. [181 114j. In the following, we concentrate 
on the part of f^s which does depend on mg. 

The part of f^g depending on rrig can be isolated through a partial derivative, which in 
turn yields a condensatel^ 

a„2fMs = (Tr[yl^])_. (2.6) 

3 \ /MS 

Thus, if we are able to measure the condensate (Tr [Ag]) lattice regularization, and convert 
the result to the MS scheme, we are able to determine the m|-dependent part of f-^ non- 
perturbatively. In principle this is, indeed, doable: the relation of the two regularization 
schemes, which is exact in the continuum limit due to the super-renormalizability of the 
theory in Eq. ()2.5p . can be found in refs. [211 [22] (see also Eq. ()3.12p below). 

Let us be a bit more specific about what we would like to achieve with the lattice simula- 
tions. Note first that in the MS scheme, the mg-dependent part of /ms is known analytically 
up to 4- loop order [23]. This corresponds to an expansion of the form 

hT[Al]) =m3 + 5ilnA + il + 4 + of4V (2.7) 
\ /MS ms m3 mg \"ig/ 

where p. is the MS scale parameter, and we have for simplicity omitted all numerical coeffi- 
cients, as well as terms containing A3. What we would like to determine is the "remainder" , i.e. 
the sum of terms beyond the level that is already known analytically. Denoting by (Tr [^o])a 
the condensate in lattice regularization, and by a the lattice spacing, the remainder is 

iTTiAl])"^ =lim|/Tr[^2A _ i _ 2 J_| _ 2 A _ ^ _ 4 , (2.8) 

\ ^ '^Vms a^olx ^ "Va a ap,) m-s 7713 ml 

where the limit inside the curly brackets takes us to the MS scheme |2H [22] , and the subse- 
quent continuum expression subtracts the known terms in Eq. (j2.7p . 

""^In the present paper we concentrate on Nc — 3, but lattice measurements have previously been carried 
out also for Nc = 2 [16]. For A^c > 4, another independent quartic coupling should be included in Eq. (12. 5|) . 

^This condensate is analogous to the Polyakov loop condensate, playing a role in various attempts at 
improved effective theories of hot QCD (see, e.g., refs. [191 120| and references therein), but it appears difficult 
to promote the analogy to a precise relation. 
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The problem with this procedure is that in practice the hmit a ^ in Eq. (j2.8p cannot 
be taken exactly, but it introduces systematic errors. To carry out the limit requires a 
fit ansatz; in three dimensions, discretization effects go like 0{a). These dominant errors 
could in principle be removed through an improvement program [21], but even if it had been 
completed (which is presently not the case) , one would still need an ansatz for the subsequent 
terms, and at this point it is for instance not clear whether logarithms like ln(l/a) should 
be included. Yet, this may have a noticeable impact on the results. Moreover, the ansatz is 
necessarily of a finite order; this means that there will remain residual 1-loop discretization 
errors in the result, while the subsequent subtraction of the continuum terms is attempting 
to take us to the 5-loop level and beyond. Evidently, this situation is unsatisfactory; for a 
demonstration of the problems encountered, see ref. |12j . 

The goal of the present paper is, then, to determine the terms corresponding to those in 
Eq. (j2.7p "exactly" in lattice regularization. That is, we do not carry out any expansion in 
ams; only one in the loop order. The corresponding result contains the counterterms needed 
for the limit in Eq. (|2.8p ; the finite terms in Eq. (|2.8p ; but also an infinite number of higher 
order corrections, starting at 0(0771,3). Then, we can write the remainder as 

(Tr [Al])^ = lim|(Tr [^^A _ ^3 ^^(^^3) _ ^2 _ ^^/^(ama) - ^fsianiA , 

\ / MS a^o L \ /a rn-g J 

(2.9) 

where fi are the functions to be determined. We of course still need to take the limit a ^ 
at the end, as indicated by Eq. (j2.9p . but we have gained in that we do not need to be as 
worried about discretization errors as before: any number yielded by the difference inside the 
curly brackets in Eq. (j2.9p is already an approximation for the remainder, and there is no 
danger of 1-loop discretization errors overtaking an interesting 5-loop continuum effect. 



3. Precise setup 

Let us now add the missing details to the formulation of the problem. The 4-loop computation 
of f-fijg in dimensional regularization has been described in detail in ref. [23j . The mass 
parameter needs renormalisation. 



m-a.bare = '^sifJ') + ^^-^ , (3.1) 

^"^^ = J^2^'^(dA + 2)(^-glX3CA + xi) , (3.2) 

where = A7rfi'^e~^^, Ca = Nc, and (Ia = — 1- Consequently the renormalised mass 
parameter satisfies the equation 

A^^^KA) = ^(dA + 2)(-5|A3Ca + Ai) . (3.3) 
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Here and in the following, we assume e — > taken in all finite quantities. It is convenient to 
define the dimensionless ratios 



y 



A3 



1 ' 



9i 



(3.4) 
(3.5) 



As indicated, we choose jl = for defining y. 

Now, it is thought that the perturbative series for the pressure of QCD converges very 
slowly [HI [7]. The reason for this can be traced back to the slow convergence of the pertur- 
bative series for 8^2 [3l [251 IS]- Making use of the results of ref. [23], and inserting the 
choice fl = g^, the known 4- loop result can be written as 



MS,At=g| 



+ 
+ 



(47r) 2 
dA 1 



+ 



Ca 



Cl 



(47r)2 4 

dA 1 
(47r)3 4yi 

+ {dA + 2 



1-2 ln(4y) 



+ x{dA + 2) \ + 



89 vr^ 11 , 

— + In 2 

12 3 3 



+ xCA{dA + 2) 



HAy) + - 



+ 



1 - ln(16y) 



-^{dA + 2f) + 



dA 1 
{AttY Ay 



1C\ 



43 491 2' 
T ~ 768^ 



+ lOxCl 



+ 



+ xCl(dA + 2 
+ 2x2CA((iA + 2) 



21n(4y) - 1 
"36-^2 



ln(4y) 



X^CAidA^lf 



TT 



x'idA + 2){dA + 8)y^ + ^'^{dA + 2f 



(3.6) 



The terms beyond Eq. (13. 6p die away as at large y, modulo possible logarithms, but 
realistic values of y are not that large (for A'c = 3, A^f = 3, y ~ 0.39 [logxo(r/Ai^) + 1.0]). 
The goal of a non-perturbative determination is therefore to sum the whole series beyond 
these known terms. 

We then move to the lattice side. To this effect, let us define the lattice action, Sa-, 
corresponding to Eq. (j2.4p . The standard Wilson discretization yields 

= /?EE{l-^I^«Tr[P,,(x)]| + 
+ 2a {tV [A2(x)] - Tr [Ao(x)C/,(x)ylo(x + i)C//(x)]} 



+ 
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+ a3^{m2,^,,TV[>l2(x)] +A3(Tr[A2(x)])'} , 



(3.7) 



where a is the lattice spacing, C/j(x) is a hnk matrix, x + i = x + ae^, where q is a unit vector, 
Pij(x) is the plaquette, and 

2Ca 



The bare mass parameter of Eq. ()3.7p reads [22j 



3, bare 



5m„ 



- 2glCA + A3(ciA + 2) 



(3i 



(3.9) 



47ra 



+^|2A3(dA + 2)(A3-5iC^A)(ln^ + C) - 25iC^A3(d^ + 2) -5 



vrS - 4(5 + p) + 2ki - K4 



(3.10) 



where S ~ 3.175911535625 is a three-dimensional hybercubic lattice integral which can be 
expressed as S = (\/3 — l)r^[^ ]r^[ii ]/487r^ [SSjfl; C,, 5, p, ki, are further lattice integrals 
which are only known numerically [22]; and fl = g^. This bare mass parameter guaran- 
tees the existence of a continuum limit for any fixed rn'^{fl) (but 0{a) discretization effects 
remain |24j ) . 

The derivative of the vacuum energy density, fa, with respect to the mass parameter (bare 
or renormalized) yields then the quadratic condensate in lattice regularization, 



(3.11) 



d^2fa = {Tr[Al] 

This can be related to the MS condensate by [HI [22] 



Tr [Al/gl] 



}im{{Tr[A'Jgi] 



levrC^ (4^)2 V 



(3.12) 



Writing the renormalised mass parameter in lattice units as 

2CAy^ 



m = amsifi = g^) 



(5rh 



/3 " 2Ca ' 

where y was defined in Eq. ()3.5p . we can express the lattice condensate as 



1 

dA 



{TT[Al/gl])^ 



1 



dAg? 



(3.13) 



(3.14) 



^We thank D. Broadhurst for bringing these references to our attention. 
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+1 



Pin \ , . . X 



+ 



[2c a) 



(t>wirh) +x4>ii{m) } + 



n=l 
3 



n=l 



(t>2nim) + 02n("l)ln 



2C2 



+ 



+ O 



2C2 



(3.15) 



Regarding the structure of this equation, we note that higher powers of logarithms than the 
terms shown do not need to be considered, as wih be explained below. 

Now, super-renormalizability guarantees that only a finite set among the functions (j)rnn^ 
4>rnn diverge in the continuum limit. In fact, as can be deduced from Eq. ()3.12p . cpoQ diverges 
as l/m and (pio diverges as ln(l/?Ti), but all the others are finite in the limit m ^ |21j . In 
this limit, they then agree with the corresponding MS scheme expressions, readily extracted 
from Eq. (j3.6p , after taking note of Eq. (j3.13|) . 

For m 7^ 0, the functions (j)mn,4'mn can be computed in lattice perturbation theory. The 
first one, (poo, follows from a 1-loop computation, while </)io,0ii require 2-loop computa- 
tions (details are given in appendix A). The functions 4)21 , 022 1 (psi^ 4'32, (pss can be de- 
duced from the fact that (Tr[Ag])a is independent of ft; on the other hand (poo, (pw, 4'ii 
have /i-dependence, emerging through the running of the MS mass parameter according to 
Eq. (13.31) . This dependence must cancel against explicit 3-loop and 4-loop logarithms, con- 
taining ln(6/a/x); these logarithms arise exclusively from the mass counterterm in Eq. (j3.10p . 
Therefore, the 3-loop and 4-loop coefficients (p2i, (P22, (psi, 4'32, (p33 can be deduced from mass 
derivatives of the 1-loop and 2-loop expressions. 

As far as the "genuine" 3-loop coefficients are concerned, we have computed explicitly only 
021 , 022 (details are given in appendix A) . These arise from graphs containing at least one 
quartic coupling, which means that most of them (with one exception) factorise into products 
of lower-order graphs. 

To display the results, we use the notation of basic lattice integrals ( Ja, la, Ha, Ga, Ba) 
explained in appendix B. Denoting furthermore 

s - 

2^- 

Ha{rh^) + G',{rn^) , (3.16) 



Kiirh') = 2 Ia{rh') + Ia{Q) 
+ 4 



/;(m2)+4(0) 



Ia{m^) + 



1 + m^drh'i 



2^ 



i:{rh^)+ia{^) 



2 + rh 9™2 



7 



+ 



+ 4 



2 + m^dri 



(47r) 



In — + C + — 



1 



2 + m-'d^2)H',{m^) + -G:im^) 



m 



+ 



laim') + ( 1 + ^ )/a(0) - — I'^im^) + 



1 



/,(m^) + - 1 + — 1,(0) 



3S 



27r- 



+ 



we obtain from appendix A, as well as from the continuum values in Eq. (|3.6p : 

1 



500 



no 



^20 
^'21 

^21 



^22 



2m 



1 



+ 0(m) , 



8irm 8tt 

m 4 4 

i(d^ + 2)[4(m2)-^]/^(rn2) 



(47r)2 

dA + 2 
(47r)2 

^2 



1 



+ 0{m) 



4 

89 11, 

ln2H hOfm) 

48 12 12 ^ ^ 



TT 



(47r)3 

{dA + 2)CAmlC3{7h^) 



{dA + 2)Ca 
(47r)3 



iln2 + ^ +0{m) 



{dA + 2)Ca . f// ~ 2\ 

2 ) 



(4^) 

(dA + 2)Ca 
(47r)3 

(dA + 2)m 



/;(m2)(ln^+C) -7i?;(m 



(4vr) 



6 



+ l-{dA + 2fm\iaim^)-^ 
Z L 47r 



2 1 

+ 2 



m 



1 6/, 



+ 



^30 



^31 



!>31 



932 



!J32 



^33 



'/>33 



dA+2 
(47r)3 



i „ln2-^^ + 0(m) 
4 16 ^ ^ 



dA + 2 
(47r)3 

(47r)4 

(4^)4 
1 

~2 



i +0{m) 



43 491 2 .^..x 
¥-1536' 



TT 



1 (dA + 2)Cl „2r^ 



(4vr) 

(47r)4 

(^A + 2)Ca 
(47r)4 



(47r 

(dA + 2)^CA ^2 

(47r)2 
(dA + 2)Ca 



1 + 0(771) 

■9 _ 7r2 _ 
4 ~ 16 ~ 



In 2 



(iA + 2 



+ 0(771) 



(rfA + 2)Ca 1 .2r- 1^2^ 
)2 2 



-m 



/^(m2)]'+[/.(m2)_^ 



(47r)4 

ciA + 2 
(47r)4 

(dA + 2)^ 

(47r)2 
(dA + 2)2 



-1 + 0(m) 



^(dA + 2)-^(<iA + 8) + 0(m) 



m 



+ 



Note that at infinite volume, /a(0) = 5]/47r, so that the functions /Ci, /C2, ^3 defined in 
Eqs. (|3.16p ~ p.l8p can be simphfied; however, in a finite volume, S/47r appearing in the mass 
counterterm is kept fixed, while /a(0), emerging from loops, gets modified (cf. appendix B). 



4. Numerical Stochastic Perturbation Theory 

In order to estimate numerically the coefficients 4'2o-, 4>3o, ^31, 4>32, 4'33j ^or which only the 
continuum values (m — > limits) are known exactly (cf. Eqs. (j3.25p . (j3.34p . (j3.35p . (j3.38p . 
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p.4ip ). we find it convenient to rewrite the action in Eq. p.7p as 
S,.u = P (l-^ReTr[P,,(x)]) - 

X, i < j 

- 2^Tr[$(x)C/i(x)$(x + i)C//(x)] + 

X, i 

+ ^|a(/5,A,m)Tr[c^2(x)]+A(Tr[cD2(x)])'| , 

where $ = y/aAo, A = oAs, m = am^{fi = g^) and, for A'c = 3, Eq. (j3.10p implies that 

^2 / 5 , \ 3.175911525625 



(4.1) 



m 



a{P,X,m) = 6<^1 + — - 6 + -A/3 



6 



lOA/3 - ^A^/?^^ In /3 + 0.08849^ + ^^^A/3 + 36.130 



87r2/32 

We write the expansion of the lattice condensate now as 



(Tr 



doo + dio^ + diiA + d2o^ + ^21 ^ + d22A^ + 



1 A A^ 

+ ^30 + dsi -^ + d32^ + dssX^ + O 

The coefficients here are related to those in Eq. (I3.15P through 



A" 



13' 



4-n 



dio 
dii 

d20 

d2i 

d22 

dso 
dsi 

d32 

dss 



2dACA4>iQ , 
dA4>ii , 

'A 



UaC^a 



m 
2dACA 
rh 

dA 
rh 
MaCI 



020 



021 + 4>2i In 



(3rh 
2CU 



'22 + <A22 In 



j3m 
2C^ 



4>30 , 



2dACA 



m 

dA 
2 



hi + In 
h2 + 4>32 In 



m 



033 + 033 In 



/3m \ 

2cuy 



(4.2) 



(4.3) 

(4.4) 

(4.5) 
(4.6) 

(4.7) 

(4.8) 

(4.9) 

(4.10) 
(4.11) 
(4.12) 
(4.13) 
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The perturbative study is concretely carried out by means of Numerical Stochastic Per- 
turbation Theory (NSPT) [271 EH]. (It would certainly also be interesting to pursue the same 
computation with standard techniques [29]; we comment on this in more detail in Sec.[6l) Its 
origins lie in Stochastic Quantization [30] , based on introducing an extra coordinate t and an 
evolution equation of the Langevin type, namely 



at$(x,i) = -5$5[$]+r?(x,t) , 



(4.14) 



where 7?(x, t) is a Gaussian noise. The usual Feynman-Gibbs path integral is recovered by 
averaging over the stochastic time, 

Z-^[ P$0[$(x)]e-^[*W] = lim - Tdt' (0[$„(x, f')])„ • (4.15) 

J t Jo ' 

In the case of gauge degrees of freedom the Langevin equation reads 

d^Ui^ = -i(vS[Ufj] +fi)Uf,, (4.16) 

where fj(x, t) is another Gaussian noise, and t is another fictitious time coordinate. Both 
time coordinates are dimensionless, being effectively measured in spatial lattice units. 

In practice, the time coordinates t, t need to be discretized as well: t = ne, t = ne, n € Z, 
with e, e ^ in the end. The discretized version of the scalar evolution, Eq. (j4.14p . reads 



<I>(x, (n + l)e) = <I>(x, ne) - e(^$(x)S' + \/er/(x, ne) , 
while the discretized version of the link evolution, Eq. (I4.16p . becomes 
Ukix, (n + l)e) = exp|-z eV^^x'S' + \/|%(x,ne) | f/fc(x. 



ne 



(4.17) 



(4.18) 



Here ri = T'^r/", fjk = T""fif.; we have rescaled the noise fields by a factor y/e, \pl ; Vfc,x = 
T"V^^, where are the generators of SU(3), normalised as Tr [T"T^] = (5"''/2; and the 
covariant derivative is defined as 

V^,x5= limi{5[e^^^"C/fc(x)] -5[C/fe(x)]} . (4.19) 
To be explicit, the expressions for the functional derivatives in Eqs. M.lTh . (I4.18P read 

+ 



5<i.(x)5 = -Y^\Uk{^)^{^^k)ul{x)^ul{x-k)^{x-k)\Jk{^-k) 



+ a$(x) + 2A$(x)Tr [<5^(x)] , 
Vfc,x5 = ^E{^fc'W-4W-|lY[p«(x)-Pi,(x) 



+ 



[c/fc(x)$(x + A:)C/,^(x),$(x) 



+ 



(4.20) 



(4.21) 
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Table 1: The masses oms = am^{j2 = g'^), counter term parts ln/3 (cf. Eq. (14. 2p ). box sizes 
N (V = a^N^), and time discretizations r studied. The box sizes were increased in unit 
steps within the intervals shown. The time step r = 5 was only used for the box sizes (a) 
N = 16, 19, 22; (b) N = 16, 19; (c) N = 16. In total, our sample consists of 298 lattices. 

Furthermore we write the gauge-field time-step in the form e = lO^^r//?, while e = lO^^r. 
To now introduce NSPT, we expand the variables as 

c^(x)^ j;5^ci>«(x), [/,(x) = l + ^r^^f (x), (4.22) 

i i=l 

where go is some small coupling; in our case, this role is played by two expansion parame- 
ters, /?~^/^ and A. This results in a hierarchical system of difference equations that can be 
numerically solved, to obtain the series in Eq. (j4.3p for (Tr [<I>^] ), for each r. Subsequently, 
we need to extrapolate to r = 0. 

Finally, we recall that the gauge field equation of motion possesses a zero-mode solution. 
When constructing the gauge field propagator, we omit this contribution; its effects are, in 
any case, insignificant in the infinite- volume limit needed for constructing Eq. ()4.3p . 

5. Data analysis 

Our approach involves three different extrapolations / interpolations in total: first, the above- 
mentioned extrapolation to r — > 0; second, an extrapolation to infinite volume (A^ oo); 
third, an interpolation between the different m = am3{fl = g"^) simulated. We discuss these 
steps one by one. The complete data sample is listed in Table [TJ 

5.1. Extrapolation r ^ 

Examples of the r — > extrapolations are shown in Fig. [H at am^ = 0.25, = 22. (We omit, 
for layout-reasons, the simplest coefficient dooO The data immediately lead to the important 
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Figure 1: The results for dy as a function of r, at am^ = 0.25, N = 22. The curves show the 
results of linear or quadratic fits. 
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observation that the shapes of the curves are practicahy the same for diQ, d2o and d^o; for 
dii, d2i and d-^i', and for ^22 and ^32. In other words, the behaviour as a function of r is 
dictated by the number of scalar coupUngs A that are associated with the coefficient. (This 
is at least partly due to the simple way in which we chose the time steps related to gauge 
and scalar field time evolutions; with some tuning, it might have been possible to optimise 
the time steps such that the time evolutions would have been more balanced [31j.) 

We have used the data at r = 5, which were the most expensive ones to produce and are 
only available for a subset of the parameter values, as a probe for the type of extrapolation 
that should be used for obtaining the r — > limits. Indeed, we choose the order of the 
polynomial fit in r low enough so that the results remain more or less stable in the inclusion 
of the r = 5 points. For dio, ^20 and dso, this requires linear extrapolations; for the other 
coefficients, we use quadratic fits. Nevertheless, we note that in some cases the results of the 
extrapolations do change by a statistically significant amount in the inclusion of the points 
with r = 5, indicating that our systematic errors may be non-negligible. 

Given the possible existence of systematic errors, it is important to crosscheck the results in 
a number of known cases. As has been discussed in Sec. El we do have exact results available, 
for any given arris and box size N, for the coefficients doo,diQ,dii,d2i,d22- In Figs. [21 [3] we 
compare the r — > extrapolations, based on 4 r's in most cases, and on 5 r's where available, 
with the exact values. 

We observe that, in general, the results of the r — > extrapolations do scatter around 
the correct values. The exception is diQ, and to a lesser extent ^21, at small volumes. We 
suspect that the reason for this is related to the way in which the zero-modes are subtracted 
in standard lattice perturbation theory (i.e. "exact values", cf. appendix B) and in NSPT, 
respectively. Nevertheless, given that the discrepancy rapidly disappears with increasing 
volume, there does not appear to be serious reason for concern. 

Moreover, we note that the extrapolations including r = 5 are in general closer to the exact 
values than those excluding it. In a few cases, the correct value is between the extrapolations 
based of 4 and 5 r's; in other words, the data point at r = 5 "overcorrects" the result of the 
extrapolation. Based on these tests, we have decided to always include extrapolations based 
on 4 and, where available, 5 r's, as independent estimates of the intercepts at r = 0. The 
extrapolations based on 5 r's have smaller error bars, and thus more weight in the subsequent 
fits; at the same time, the inclusion of the extrapolations based on 4 r's allows us to correct 
for the mentioned overshooting in the cases where it does take place. 

5.2. Extrapolation N ^ 00 

Given the results of the r ^ extrapolations (the complete data set, apart from doQ, is shown 
in Fig. Hj), the next step is to extrapolate to infinite volume. 

As Fig. m shows, finite-volume effects become small at large volumes. However, the box 
size N = L/a required for this grows as the mass am^, decreases (the behaviour is more or 
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Figure 2: The results for doo and dio after the r — > extrapolation, as a function of the 
inverse box size in "physical" units {m^L = am^N), together with the exact values. The 
(infinite volume) "extrapolation" denoted with the closed square refers to the procedure 
defined around Eq. (15. 3p . 



less universal as a function of m^L). In addition, some of the coefficients appear to require 
larger volumes than others. For the smallest masses, in particular, we are in many cases not 
yet in a region where all volume dependence has died out. 

To be able to deal with this situation, some theoretical knowledge about the functional 
dependence on the finite volume is needed. The situation is complicated by the fact that there 
are both massless and massive fields in the system; therefore both powerlike and exponential 
volume dependences appear. However, an inspection of the known cases (Figs. O [3]) suggests 
than in practice the magnitude of the power corrections is much smaller than that of the 
exponential ones so that, strangely enough, the latter dominate in the volumes where our 
data lies. In this situation, we could then expect the dominant volume behaviour to be some 
exponential, dij{ra^L) — dij [oo) ~ exp(— m3L)/(m3L)". Unfortunately, an inspection of some 
of the known cases (particularly ^22) shows that, again because of the fairly small volumes 
reached in practice, the behaviour is not given by a simple exponential, but that there are 
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Figure 3: Like Fig. [2] but for du, d2i, ^22- 
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at least two competing universal functions, because the position of the maximum in m^L 
evolves slightly with arris- At the same, allowing for too many free functional forms in the 
fits does not lead to good results either, because a long extrapolation is needed, so that the 
ansatz needs to be fairly constrained. 

Having tested many procedures in the cases where exact results are available, we have 
finally chosen the following strategy in order to deal with these challenges. Let us consider 
a mass like am^ = 0.8. Fig. S] shows that this one has a reasonable plateau at affordable 

< 15 for all the coefficients dij, but still the behaviour of the data is not too flat (i.e., some 
volume dependence is detectable). One can then extract an infinite- volume value dij (oo) by 
fitting a constant to data in the range of the plateau, and subtract it from the data in order 
to obtain the quantities 



Subsequently, one can try to obtain a reasonable interpolating fit fij{msL) for gij{rasL), 
allowing to go also to other values of m^L than those simulated at am^ = 0.8. In practice, 
we find that in the range m^L > 2.5 that we have considered (cf. Table [T]), our data for 
am3 = 0.8 can be well modelled, for instance, by the ansatz 



where are fit parameters. We have also experimented with other ansatze, but do not 
find a significant effect on our final results. 

After this empirical determination of the finite-volume effects in one well-controlled case, 
we can go back to the other masses am'^, and use the fits fij as a constrained ansatz. However, 
as already mentioned, there are cases, such as ^22, where the position of the maximum evolves 
with ams; therefore the results cannot be described by one universal function in our modest 
volumes. To incorporate this fact, we allow Eq. (15. 2p to in general split up into two functions, 
and take a finite-size scaling ansatz of the form 



where x' = m'^L; dij{x') are the direct measurements at the mass am'^ for various A^ = L/a] 
and dij(oo), Aij{am'^) and Bij{am'^) are volume-independent fit coefficients. Among the 
exactly known coefficients, the only case where the results change significantly while going 
from Eq. (j5.2p to the more general Eq. (j5.3p is precisely ^22 (00); among the unknown ones, 
the ansatz does systematically affect also ^20 (00), and particularly (£30(00), at the smallest 
masses. For instance, in the last case, the values of ^30(00) would be as much as ~ lOo" higher 
at the two smallest masses if we employed Eq. (j5.2p throughout. In the following, we cite 
results based on Eq. (j5.3p . for reasons now to be explained. 

In Figs. El [31 the results of such fits are compared with the exact results at the small- 
est masses. We do find compatibility within statistical errors (~ 2a) in all cases. The 



gij{msL) = dij{msL) - dij{oo) . 







(5.3) 
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results for the coefficients dij = dij{oo), both fitted and exact, are given in Table [2] for all 
masses. Only one "failure" can be detected, namely the coefficient dio at the largest masses 
am^ = 0.80, 1.00, where our (very small) NSPT error bars appear to be underestimated (the 
difference is ~ lOa). Given that the largest masses are the least important ones for the 
subsequent steps, we have decided to let this problem "pass"; let us stress, in any case, that 
the overall excellent agreement is a very non-trivial result, as a long extrapolation needs to 
be carried out, and could only be achieved after a considerable amount of experimenting with 
various procedures. 

Encouraged by these tests, as well as by the indication in Figs. [H U that d2o, d^Q could 
more or less behave like dio; d^i like d2i and dn; and (^32 like ^22; we then apply the same 
procedure to the remaining coefficients. The results are given in the lower-most panel in 
Table [2] (errors are statistical only). Finally, all the results, but with the normalization of 
Sec. m are shown in Fig. [5] (the exactly known values of (j)2i, 4>22-, 4>?,i-, 4'32, 033 have been 
used as input to convert ^21, ^22, ^31, c?32, c?33 to (j)2i, 4>22, (psi, 032, 4>33, respectively; cf. 
Eqs. (j4.8p - (j4.13p . The functions 021, 022 , 03i, 032, 033 have been numerically crosschecked 
only at am^ = 1.00, where simulations with two different ln/3 were carried out; cf. Table [1]). 



5.3. Interpolation in am^ 

The remaining task is to provide interpolating fits for our functions (pij (arris), 0ij (0^77-3). 
Indeed, lattice simulations such as those in ref. [12] correspond to values of am^, given by 
Eq. (j3.13p . as well as lattice spacings In/?, which are in the range of our study, but seldom 
coincide exactly with our values. The purpose of the interpolating fits is to nevertheless make 
our results usable for the analysis of lattice simulations. 

In order to carry out the interpolating fits, a fit ansatz is again needed. Since the continuum 
limit corresponds to am^ ^ 0, it may be reasonable to use a finite-order polynomial in am^ 
for this purpose. However, individual graphs do lead to other structures as well, particularly 
logarithms like 07713 ln(l/am3) (cf. Eqs. ()B.9p . ()B.12p ). Even though these logarithms cancel 
in all the analytically known terms, we are not aware of a proof excluding them in general. In 
particular, the lattice simulations of ref. [12j strongly suggest the presence of such a logarithm, 
affecting the approach to the continuum limit, and it would then be natural for it to appear 
in the coefficient 02o, which is numerically the most important unknown ingredient entering 
the analysis of ref. |12j . 

Given these considerations, we have carried out fits of two types. Defining 



^00 



^10 



87ram3 
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Figure 5: Shown are the coefficients in Eq. (j3.15p . as a function of am^. The polynomial fits 
are of third order in ams, and have been constrained to go through the continuum points; 
the logarithmic fits include the additional term am^\n.{\ / amz) , cf. Eq. (|5.8p . 
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we have considered a polynomial fit, 

(j)ij{am3) = a\f + a^f x arris + a^f x {amsf + af^^ x {amsf , (5.7) 
and a logarithmic one, 

(pljiam^) = + hfP X ams In ^ + h^^ x ams + ftjf x {am^f + b\f x {am^if . (5.8) 

The fits have been constrained to have the correct continuum values aij{0), bij(0) which are 
known in all cases (cf. Sec. [3]). The fit functions are illustrated in Fig. \5\ and the results for 
the coefficients are given in Tables [3l HI 

Two main observations can be made from the fits: 

• In all cases where exact results are available, the logarithmic fits agree reasonably well 
with the polynomial ones. This is expected in the sense that we know that no logarithms 
exist in the exactly known functions. 

• For the most important unknown function, <j)2o, the logarithmic fit does appear to pro- 
duce a markedly better description of our data than a polynomial fit; x^/d.o.f. decreases 
dramatically, from ~ 55 to ~ 0.13. (This is the case also for the second-most impor- 
tant unknown, 03o, where x^/d.o.f. decreases from ~ 20 to 1.2.) This observation 
would appear to be in accordance with the indications from lattice simulations |12j . 
In fact, the authors in ref. [12] estimated the logarithmic term to be (Tr[^Q/(7|]) ~ 

... + (0.10... 0.13) X5r|aln(l/a), which in our units converts to b'^^p - (0.10 ... 0.13)/^^ = 

(1') 

0.012 ... 0.017. Though the agreement with our value, 620 ~ 0.00973, is not perfect, 
the order of magnitude is the same. It would be interesting to re-analyze the results of 
ref. [12] with the coefficient of the logarithm fixed to our value. 

In summary, then, there appear to be good reasons to expect the presence of logarithms in (j)2o 
and (psQ. To unambiguously confirm this expectation, it would obviously be very interesting to 
find a way to improve on the accuracy of the determination of these coefficients, controlling 
in particular the difficult-to-estimate systematic errors that are related to the r ^ and 
— > 00 extrapolations in the present method. It would also be important to improve on 
the values of the coefficients associated with scalar self-couplings, (/>32, (^'33, although 
this would mostly serve as a theoretical consistency check, given that the values of Xs/g^ 
corresponding to physical finite-temperature QCD are very small |16j . 

6. Conclusions 

The purpose of this paper has been to estimate the (Debye) mass dependent part of the 
vacuum energy density of the three-dimensional SU(3) + adjoint Higgs theory, up to 4- 
loop order in lattice perturbation theory. The result can be parametrized in terms of the 
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coefficients 'Pij^ defined in Eq. (j3.15p . We have worked out the expressions for a number 
of these coefficients analytically (Eqs. (|3.19p - (j3.42p ). and estimated the remaining ones, 4>2o, 
<;^30; <^3i; </'32; 0337 for which only the continuum values are known analytically, with the 
help of Numerical Stochastic Perturbation Theory. The results are illustrated in Fig. [5l and 
parametrized in terms of simple fits in Eqs. (15.71) . (15. Sp . 

The main practical use of our results is that when combined with lattice Monte Carlo data, 
they should allow to improve on the analysis of the sum (beyond the known 4-loop order) 
of infrared sensitive contributions to the pressure [25] and quark number susceptibilities [12] 
of hot QCD, given that discretization errors up to the 4-loop order can now be subtracted. 
However, our results might also have some theoretical interest beyond these particular appli- 
cations. For instance, they serve as a consistency check of the 4-loop MS-scheme computation 
of ref. [23] (in the sense that our results appear to be consistent with the continuum values 
indicated in Fig. [5]), as well as of the super-renormalizability of the theory considered and 
the power-counting arguments presented for it in ref. [21] (in the sense that no indications of 
ultraviolet divergences apart from the known 1-loop and 2-loop ones in 0oo! 9^10 were seen). 

Concerning the new coefficients (/)2o, </'3o> </'3i; <^32 and (/>33, we unfortunately have to ac- 
knowledge that it appears difficult to improve significantly on the accuracy of our results with 
the present techniques. The problem is that two different extrapolations, r ^ and N ^ oo, 
are needed in order to obtain the values at any fixed ama, and both of these extrapolations 
introduce systematic and statistical errors. Therefore there is a need to crosscheck our results, 
and improve upon them, with standard techniques [29]. Nevertheless, even in that approach, 
our study should serve as a basic framework. In particular, we would like to stress the insight 
that it is valuable to determine the coefficients <j)ij, (pij as functions of am^, rather than to 
carry out an expansion in small arris, since realistic values of ams (~ 0.1 .. . 0.5 |12j ) are in 
a region where the functions show more structure than just linear terms (cf. Fig. [5]). The 
most important coefficients to determine are (j)2o and cp^o, which are independent of the scalar 
self-coupling, and for am^ ^ 1, a concrete challenge is to confirm or disprove the existence 
of the logarithmic term ~ 0{alii{l/a)) in (f>2o, for which independent indications have been 
seen in ref. [12] and in the present study. 
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Appendix A. Graph-by-graph results for the vacuum energy density 



We list here results for the mass-dependent part of the vacuum energy density fa of the 
theory in Eq. (13.7p (through Eq. ()3.1ip this produces the condensate that we are interested 
in). At the 3-loop level we have only kept terms involving at least one A3. In the graphical 
notation to be used below, solid/wiggly lines represent tree- level Aq/ Ai lattice propagators, 
respectively. 



1-loop and 2-loop graphs. For brevity, we denote in the following m = ms^fi). The 
integrals appearing are defined in appendix B, and in terms of these, the results read (d = 3): 



O - T^"(- 



CO 



A.. 



(A.l) 

'-dA{dA + 2)[Ia{m')]\ (A.2) 
© + CO " ^dACA{{2d - 4)/,(0)/„(m2) + [Ia(m2)]2 + Am'Hair 



3-loop graphs involving A3 

ceo 

(DO-CCO 



-fdAidA + 2f[Ia{m^)flUm') 



X2 
-^3 



dA{dA + 2)Ba{m^ 



(A.3) 

(A.4) 
(A.5) 



53A3 



dACA{dA + 2)Ia{m^){{2d-4)Ia{0)l'aim^) + (A.6) 



+ 2Ia{m'K{m') + AHa{m') + Am'H'^{m') + 



. 



4(0)4(m2) + m'laW^im') + G',{m')] } 



(A.7) 



Mass counterterm contributions up to 3-loop order. Gauge field "mass countert- 
erms" (arising from the Haar integration measure) are not displayed, as they do not contribute 
to terms involving at least one A3. 

dA 



[m 



A3 



00 



dA{dA + 2)6mllaim')l',{m^) , 



^ ^-}dACASml{i2d - A)UO)I',im') + 



(A.8) 
(A.9) 
(A.IO) 
(A.ll) 
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+ a" [/,(0)4(m2) + m24(0)/;(m2) + G',{m')\ } 



Appendix B. Basic lattice integrals 



We detail here the definitions of the basic lattice integrals that appear in the expressions 
discussed in Sec. [3l The integration measure is 

'-^ d^p 



dp 



. (2vr)3 



and we define the standard lattice momenta as 



2 sm — 
2 



i=l 



(B.l) 



(B.2) 



The integrals appearing then read (we use hats as a reminder of the use of lattice units in 
these expressions): 



dp \n(jP' + m 

_|_ 

dpdg ■ 



dpdq 



{p^ + ■im?){(p + ■im?){p + q) 



2 ' 



(B.3) 
(B.4) 
(B.5) 

(B.6) 



(p2 + rn?){(p + rh?){p + 

Small-m expansions for these functions have been worked out in refs. |2H [22] and are given, 
to the order that was used for the small-m expansions in Sec. [31 by 

J^(m2) = l{rn^)^ (B.7) 

T.-m + 0{'m^) , (B.8) 



1 



(47r)2 
1 



3 1 

In — + - +C + 0{rh) 
m 2 

16K1 -4(5m2 + 0(r7^3) 



(B.9) 
(B.IO) 



(47r)2 L 

where S, C, ki and 5 are numerical coefficients mentioned below Eq. (|3.10|) . Furthermore, we 
define a 3-loop "basketball" integral through 



dp dq dr ■ 



1 



(B.ll) 
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In this case the small-m expansion reads 

1 



S X 6*- m( 41n — + 4C + 6 ) + Ofm^) 
2m 



(B.12) 



The derivation of this result, which has to our knowledge not appeared in the literature 
before, as well as a numerical estimate for the new coefficient 9, are given in appendix C. 
The expressions in appendix A employ the functions 







(B.13) 




a 

= Hairh^), 


(B.14) 
(B.15) 
(B.16) 




a 


(B.17) 



where a is the lattice spacing, and rh = am. 

In a finite volume, V = (aN)^, the momentum integrations get replaced with 



N-l 7V-1 



Ap/(p) - ]^ E E E /(^) . 



(B.18) 



ni=0 n2=0 n3=0 

where n = (ni,n2,n3). In the case of massless propagators, the zero-mode is left out. 



Appendix C. The 3-loop basketball in lattice regularization 



In order to work out the expansion in Eq. ()B.12p . we find it convenient to return from lattice 
units to physical units, considering then the function in Eq. (|B.17|) . Let us introduce the 
scalar propagator on the lattice. 



Da{x;m 

where now p'^ = Yli=iPh Pi ~ 
continuum limit a — > 0, 



(C.l) 



DQ{x]m) 



For X 7^ 0, the propagator remains finite in the 
exp(— m|a;|) 



47r|x| 



while for x = 0, it contains a linear divergence as shown in Eq. (jB.Sp . 

1 



Da{Q]m) 



E — m + Oirh) 



(C.2) 



(C.3) 
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The integral we will be concerned with here is of the "basketball" type, 

n 



(C.4) 



For n = 3, this equals the integral Ha defined in Eq. (|B.5p . and for general masses, the 
expansion close to the continuum limit is of the form [21] 



1 



(4vr) 



6 1 
In^^ + - +C + 0{ami) 



(C.5) 



The challenge now is to find the corresponding expansion for n = 4. 

The basic approach we follow is similar to the one used for the basketball integral in 
dimensional regularization in refs. [32]. The summation over configuration space in Eq. ()C.4p 
is divided into two regions, |x| < r and |x| > r. We assume ami 1, and can thus choose 



1 

a <C r < — . 



(C.6) 



In the region |x| < r, we now have \x\mi <^ 1, and can expand in the masses; in the region 
|x| > r, we have a/\x\ <^ 1, and can use the continuum approximation for the propagators. 

The region \x\ > r. The region of large |x| gives a contribution which remains finite in the 
limit a ^ 0. It can thus be evaluated employing Eq. ()C.2p . The integral is elementary, and 
we obtain 



limAu|>,i?(4)({ma) 



a->0 



47r|xpd|x| JJl)o(a;; 



rrii 



i=l 



(47r)= 



- + M(lnMr + - l) 



+ 0{M^r) , (C.7) 



where M = mi + m2 + + 7714. Note in passing, for future reference, that 

3 



47r|xpd|x| JjL>o(a;; 



1 



rrii 



i=l 



In Mr - 7e + 0{Mr) , 



(C.J 



where M = mi + m2 + m^. 



The region |x| < r. In the region of small we rewrite the propagator of Eq. (jC.ip in 
the equivalent form 



Da(x;m) 



n/a ^; 



P 



1 



1 



m 



p2 p2 _|_ ^2 p2 p2 (^p2 _|_ j^2 



The point of this rewriting is that for a ^ and x 7^ 0, the first term behaves as ~ l/47r|x|, 
the next two combine into a constant ~ —m/in, while the last term, which is both ultraviolet 
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and infrared finite, behaves as 0(m^|x|). Therefore, if we consider the hmit a ^ 0, rnii <C 1, 



the last term does not contribute in AM^^Bji'\{mi}): 



The same is true for the case that there are two or more appearances of the constant term: 

\xfd\x\ .^ f ..^ 0{mi)0{mj) ~ 0(mim,r) . (C.ll) 

Thus the only contributions come from four appearances of the first term in Eq. ()C.9p . and 
three appearances of the first term as well as one appearance of the middle terms: 

limA|,|<,i?(^)({mJ) = limTa^[Da{x;0)f (C.12) 

a— >0 a^O ^ — ' 

\x\<r 

+ lim|j^[z),(0;mi) -Z)a(0,0)] ^ a^Daix;0)A + Oim^r) . 

^i=l \x\<r ^ 

Here and in the following, lima^o is meant in a symbolic sense, since the sums actually 
diverge as a — > 0. 

Given that Da{0; mi) — Da{0, 0) is known (cf. Eq. ()C.3P ). we are left with the evaluation of 
the sums on the first and second rows in Eq. ()C.12p . Since the propagators are massless, the 
outcomes only depend on the ratio r/a, where a <^ r. The first of the sums can be performed 
by extending the sum to be over all space, and taking the continuum limit in the resulting 
subtraction, which is ultraviolet finite: 

hm V [Da{x; 0)]^ = lim V [Da{x; 0)f - lim V [Z)„(x; 0)]^ 

a— >0 ^ — ' a^O ^ — ' a^O ^ — ' 

\x\<r X \x\>r 

= limVa3[Z),(x;0)]4--i^-, (C.13) 

where we usedEq. (lOTl) . The latter sum is slightly more difficult because it would be infrared 
divergent at large distances, but it can be performed as above, once we regulate the infrared 
with mass terms, and use then Eqs. (jC.Sp . (|C.8P : 

,3 3 . 

lim V [Da{x] 0)f = lim <^ lim V TJ Da{x; rm) - lim V TJ Daix; m,) \ 

|x|<r ' ^ X i=l \x\>r i=l ' 

= limi— i-^(ln-^ + i +C + lnMr + 7E ) I . (C.14) 
a^o\(4^)2 \ aM 2 ^ ' J I 

The infrared regulator M is seen to cancel in Eq. ()C.14p . as it should. 
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Combining now Eqs. ([03]) . (ja7|) . ([012]) . (jCTTS]) and ([CHl) . terms singular in r cancel, 
and the limit rriir can be taken. The Euler gamma-constants 7e also cancel, and we 
finally obtain 

Bi'Him,}) = Y,a' [Daix; 0)f - ^ (^In ^ + ^ + c) + 0(am?) . (C.15) 

The task remains to evaluate the sum in Eq. (|C.15p . Employing the techniques introduced 
in ref. [S^ and worked out for the three-dimensional case in ref. [M], we find that 

Y^a^D,ix;0)]' ^ -j^e , = 3.0122(1), (C.16) 

where the number in parentheses indicates the uncertainty of the last digit. Combining 
Eqs. (|UT5]) . (|UT6]) leads to Eq. (|Bl^ . 

Finally, recall for completeness that in dimensional regularization (DR) at (i = 3 — 2e, the 
2-loop [35J and 3- loop [3B] basketball integrals read (M = Yli=i mi, M = X]^=i ^i) 

For a complete discussion of basketball integrals in dimensional regularization, see ref. [37] . 
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-2.16(6) 


-0.060(19) 


0.30 


14.09(3) 


81.8(3) 


9.6(3) 


-1.73(5) 


-0.038(10) 


0.40 


10.87(1) 


52.2(1) 


9.4(1) 


-1.27(1) 


-0.029(3) 


0.50 


8.74(1) 


36.7(1) 


8.61(5) 


-0.916(7) 


-0.0201(7) 


0.60 


7.18(1) 


27.2(1) 


7.59(5) 


-0.670(6) 


-0.0155(5) 


0.80 


5.058(2) 


16.11(2) 


6.035(7) ■ 


-0.3522(9) 


-0.01423(6) 


1.00 


3.664(1) 


10.16(1) 


4.652(4) ■ 


-0.1605(5) 


-0.01327(3) 


1.00* 


3.666(6) 


10.19(6) 


6.46(2) 


-0.196(3) 


-0.01677(11) 



Table 2: The coefficients dij (cf. Eq. ()4.3p ) in the infinite- volume limit. The starred mass 
refers to In /3 = In 80; the others to In /? = In 24. The numbers in parentheses indicate the 
statistical errors of the last digits shown. 
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coefRcients 










xVd.o.f. 


*i 




-0 039789 


-0 0062086 


0056620 

yj*yjyjKjyjyj^yj 


-0 00076693 




*i 




038310 


-0 014881 


-0 00082279 

yj *yjyjyjK^ ^ ^ 1 


0029967 




*i 


hi 


015831 


0076464 

yj ryjyj 1 yj^yj^ 


-0 0094444 


0017264 




C 


P20 


0092578 


027798 

yj myj ^ 1 1 


-0 039072 


014739 


55.07 




hi 


0.010909 


0.014947 


-0.0022681 


-0.0043976 






t>22 


-0.0053827 


-0.0050144 


0.0074469 


-0.0014995 






ho 


0.0024038 


0.0071490 


-0.0020206 


-0.0016577 


19.83 


( 


hi 


0.0013885 


-0.0035354 


0.024195 


-0.013013 


0.996 




h2 


-0.0018767 


-0.0049224 


0.0013548 


0.0029473 


1.966 




hs 


-0.00031675 


-0.00036476 


-0.00032547 


-0.00014179 


1.063 




hi 


0.0075590 


0.0023970 


-0.0033645 


0.00069377 






h2 


-0.0025197 


-0.00079883 


0.0011212 


-0.00023111 






hi 


0.0036091 


0.0031005 


-0.00069559 


-0.00088573 




*<-. 


h2 


-0.0012030 


-0.0026312 


0.0043532 


-0.0011224 






hs 


0.0 


0.00053248 


-0.0013736 


0.00047246 





Table 3: The fit coefficients allowing to estimate the functions (pij, (pij in the range 0.0 < 
ams < 1.0, according to Eq. (j5.7p . In the cases marked with a star, the fits have been carried 
out to the exact results rather than to NSPT data (we cite no x^/d-o.f. here because no 
error bars can be assigned to the exact numbers); the accuracy of their description through 
a polynomial fit is on the per cent level. 



coefficients 






b^^ 




6^^) 


xVd.o.f. 


020 


0.0092578 


0.0097310 


0.0090404 


-0.010317 


0.0047392 


0.128 


030 


0.0024038 


-0.0036880 


0.014462 


-0.013444 


0.0024569 


1.185 


031 


0.0013885 


-0.0017147 


-0.00022591 


0.019122 


-0.011248 


0.932 


032 


-0.0018767 


0.0058403 


-0.016324 


0.018943 


-0.0032415 


0.361 


033 


-0.00031675 


-0.0029965 


0.0052048 


-0.0086842 


0.0026476 


0.737 



Table 4: The fit coefficients allowing to estimate the functions cj)20i 030; 03i> 032) (Ass in the 
range 0.0 < am^ < 1.0, according to Eq. ()5.8p . 
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